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We analyze the performance of the Resonator/zero- Qubit (RezQu) architecture in which the 
qubits are complemented with memory resonators and coupled via a resonator bus. Separating 
the stored information from the rest of the processing circuit by at least two coupling steps and 
the zero qubit state results in a significant increase of the ON/OFF ratio and the reduction of the 
idling error. Assuming no decoherence, we calculate such idling error, as well as the errors for the 
MOVE operation and tunneling measurement, and show that the RezQu architecture can provide 
high fidelity performance required for medium-scale quantum information processing. 
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I. INTRODUCTION 

Superconducting circuits with Josephson junctions are 
steadily gaining attention as promising candidates for the 
realization of a quantum computer [1]. Over the last sev- 
eral years, significant progress has been made in prepar- 
ing, controlling, and measuring the macroscopic quantum 
states of such circuits [2-12]. However, the two major 
roadblocks - scalability and decoherence - still remain, 
impeding the development of a workable prototype. The 
Resonator/zero- Qubit (RezQu) protocol presented here 
aims to address these limitations at a low-level (hard- 
ware cell) architecture [13]. 

A RezQu device consists of a set of superconducting 
qubits (e.g., phase qubits [14]), each of which is capaci- 
tively coupled to its own memory resonator and also ca- 
pacitively coupled to a common resonator bus, as shown 
in Fig. 1 [7, 15-17]. The bus is used for coupling oper- 
ations between qubits, while the memory resonators are 
used for information storage when the logic qubits are 
idling. With coupling capacitors being fixed and rela- 
tively small, qubit coupling is adjusted by varying the 
qubit frequency, which is brought in and out of reso- 
nance with the two resonators. For a one-qubit opera- 
tion, quantum information is moved from the memory to 
the qubit, where a microwave pulse is applied. A natural 
two-qubit operation is the controllcd-Z gate, for which 
one qubit state is moved to the bus, while the other qubit 
frequency is tuned close to resonance with the bus for a 
precise duration [6, 18-20]. 

Most importantly, the information stored in resonators 
is separated from the rest of the processing circuit by the 
known qubit state |0) and at least two coupling steps, 
thus reducing crosstalk error during idling. Also, the 
problem of spectral crowding is essentially eliminated be- 
cause the two-step resonance between empty qubits is 
not harmful, while the four-step coupling between mem- 
ory resonators is negligible. Therefore the resonator fre- 
quencies, which are set by fabrication, can be close to 
each other, decreasing sensitivity to phase errors in the 
clock. Thus the RezQu architecture essentially solves the 
inherent ON/OFF ratio problem of the fixed capacitive 
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FIG. 1: Schematic diagram of the RezQu architecture: m - 
memory resonators, q - qubits, b - bus. We assume frequen- 
cies ~7 GHz for the memories, ~6 GHz for the bus, and qubit 
frequencies are varied between these values. 



coupling without using a more complicated scheme of a 
tunable coupling [8, 21, 22]. As an additional benefit, in- 
formation storage in resonators increases coherence time 
compared to storage in the qubits. We note that the 
idea of using resonators to couple qubits has been sug- 
gested by many authors [23-32]. The use of resonators 
as quantum memories has also been previously proposed 
[31, 33-35]. However, putting the two ideas together in a 
single architecture results in the new qualitative advan- 
tages, which have not been discussed. 

In this paper we briefly consider the relation between 
the logical and the physical qubit states and then ana- 
lyze several basic operations in the RezQu architecture. 
In particular, for a truncated three-component memory- 
qubit-bus RezQu device we focus on the idling error, in- 
formation transfer (MOVE) between the qubit and its 
memory, and the tunneling measurement. The analysis 
of the controlled- Z gate will be presented elsewhere [36]. 
For simplicity decoherence is neglected. 



2 



II. LOGICAL VS PHYSICAL QUBITS 

We begin by recalling an important difference between 
logical and actual, physical qubits. The difference stems 
from the fact that in the language of quantum circuit di- 
agrams the idling qubits are always presumed to be sta- 
tionary [37], while superconducting qubits evolve with a 
6-7 GHz frequency even in idling, which leads to accu- 
mulation of relative phases. Also, the always-on coupling 
leads to fast small- amplitude oscillations of the "bare"- 
state populations during the off-resonant idling. A nat- 
ural way to avoid the latter problem is to define logical 
multi-qubit states to be the (dressed) eigenstates of the 
whole system (see e.g., [5, 26, 28, 29, 38, 39] for discussion 
of the dressed states). Then the only evolution in idling 
is the phase accumulation for each logic state. However, 
there are 2 N logic states for N qubits, and using 2 N "lo- 
cal clocks" (rotating frames) would require an exponen- 
tial overhead to calibrate the phases. The present-day 
experiments with two or three qubits often use this un- 
scalable way, but it will not work for TV > 10. A scalable 
way is to choose only N rotating frames, which corre- 
spond to one-qubit logical states, and treat the frequen- 
cies of multi-qubit states only approximately, as sums 
of the corresponding single-qubit frequencies. The use 
of such non-exact rotating frames for multi-qubit states 
leads to what we call an idling error, which is analyzed 
in the next section. 

Notice that it is sufficient to establish a correspondence 
between logical and physical states only at some time 
moments between the gates. Moreover, this correspon- 
dence may be different at different moments. At those 
moments, the bus is empty, the 27V + 1 system compo- 
nents (N qubits, N memories, and the bus) are well- 
detuned from their neighbors, and for each logical qubit 
it is unambiguously known whether the corresponding 
quantum information is located in the memory or in the 
qubit. Thus the eigenstates corresponding to the logical 
states are well-defined and the physical-to-logical corre- 
spondence is naturally established by projecting onto the 
2 N computational eigenstates, while occupations of other 
eigenstates should be regarded as an error. In the sim- 
plest modular construction of an algorithm we should not 
attempt to correct the error of a given gate by the sub- 
sequent gates. Then for the overall error we only need to 
characterize the errors of individual gates, as well as the 
idling error. 

One may think that defining logical states via the 
eigenstates of the whole system may present a techni- 
cal problem in an algorithm design. However, this is not 
really a problem for the following reasons. First, we need 
the conversion into the basis of eigenstates only at the 
start and end of a quantum gate, while the design of a 
gate is modular and can be done using any convenient 
basis. Second, in practice, we can truncate the system 
to calculate the eigenstates approximately, making sure 
that the error due to truncation is sufficiently small. Sim- 
ilar truncation with a limited error is needed in practical 



gate design. 

As mentioned above, the physical-to-logical correspon- 
dence rule can be different at every point between the 
gates. For the correspondence based on eigenstates we 
are free to choose N single-excitation phases arbitrar- 
ily. In spite of this freedom, for dcfinitcness, it makes 
sense to relate all the single-excitation phases to a par- 
ticular fixed time moment in an algorithm. Then a shift 
of the gate start time leads to easily calculable phase 
shifts, which accumulate with the frequencies equal to 
the change of single-excitation frequencies before and af- 
ter the gate. Such shift is useful for the adjustment of 
relative single-excitation phases [40] . Another way of the 
single-excitation phase adjustment is by using "qubit fre- 
quency excursions" [7]. The ease of these adjustments 
significantly simplifies design of quantum gates, because 
we essentially should not wary about the single-excitation 
phases. 

Notice that initial generation of high-fidelity single- 
excitation eigenstates is much easier experimentally than 
generation of the bare states. This is because a typi- 
cal duration of the qubit excitation pulses is significantly 
longer than inverse detuning between the qubit and res- 
onators. We have checked this advantage of using eigen- 
states versus bare states numerically in a simple model 
with typical parameters [7] of a RezQu device (the er- 
ror decrease is about two orders of magnitude). Simi- 
larly, the standard one-qubit operations essentially oper- 
ate with the eigenstates rather than with the bare states. 

III. IDLING ERROR 

Before discussing the idling error in the RezQu archi- 
tecture, let us consider a simpler case of two directly cou- 
pled qubits. Then in idling the wavefunction evolves as 
\ip(t)) = aoo_e-" oot |00) +a ie- 4£oit |0T) +a 10 e- 4£lot |T0) + 
ane~ 4£ll *|ll), where we denote the logical (eigen)states 
with an overline, their corresponding (eigen)energies with 
ey, and amplitudes at t = with a^. However, for the 
desired evolution | ^desired (*)) the last term should be re- 
placed with a 11 e _i ( eol+eio_eo °)*|TT); then only two rotat- 
ing frames (clocks) with frequencies eoi — eoo and eio — £oo 
are needed. We see that the phase difference accumulates 
with the frequency Qzz = (en — £oi) — (^io — £oo)> an d 
therefore the idling error due to qubit coupling accumu- 
lates over a time t as 

Err = l-|(V>de S ired(*)|V>W>| 2 

^ Ki| 2 (i - Ki| 2 )0W) 2 < (Clzzt) 2 , (1) 

where we assumed Clzzt <C 1. (The frequency Qzz is 
defined in the same way as in Ref. [39] for a two-qubit 

cr^V^ interaction.) The error is state-dependent, but 
in this paper we will always consider error estimates for 
the worst-case scenario. 

In the RezQu architecture, the main contribution to 
the idling error comes from interaction between a mem- 
ory resonator, in which quantum information is stored, 
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and the bus, which is constantly used for quantum gates 
between other qubits. By analogy with the above case, 
for the truncated memory-qubit-bus (mqb) system the 
idling error can be estimated as 

Err ~ (ttzz t) 2 , fizz = (eioi - eooi) - (eioo - eooo), (2) 

where the eigenenergies correspond to the logical eigen- 
states 1 101), |00T), |T00), |000), and in our \mqb) nota- 
tion the sequence of symbols represents the states of the 
memory resonator, the qubit, and the bus. Notice that 
the qubit here is always in state |0), and fizz is essen- 
tially the difference between the effective frequencies of 
the memory resonator in the presence and absence of the 
bus excitation. 

To find fizz we use the rotating wave approximation 
(RWA) ; then the dynamics of the mqb system is described 
by the Hamiltonian (we use H = 1) 



H(t) = 





u q (t) 

2u q (t) - rj 



+ w m a4a m + uj b a\a b 



+ 



(a} n a q + a m cr+) + g b (a q a\ + a+a^j , 

+ 9d (ali a b + a m alj , (3) 

where the qubit frequency u q may vary in time, while the 
qubit anharmonicity 77 is assumed to be constant, 
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(4) 



ujb are 



are the qubit lowering and raising operators, 
the memory and the bus frequencies (which are presumed 
to be fixed), a^, a m , at, a b are the creation/annihilation 
operators for the memory and the bus photons, and g m , 
g b are the memory-qubit and qubit-bus coupling con- 
stants. The last term in Eq. (3) describes the direct 
(electrostatic) memory-bus coupling; replacing a qubit 
in Fig. 1 with a lumped tank circuit it is found [41] to be 



9d = 2g m gb/u q . 



(5) 



It is typically smaller than the effective memory-bus cou- 
pling via the virtual excitation of the qubit because the 
detunings \u m — uj q \ and \tu b — ui q \ between the elements 
are much smaller than their frequencies; because of that, 

we often neglect g^. From the physical model it is easy to 

1/2 

show that g m and g b are proportional to u> q and there- 
fore change when the qubit frequency is varied; however 
for simplicity we will assume constant g m and g b . 

Neglecting g^, in fourth order we find (see Appendix), 



fizz 



-^9 1 m9 b r t u m + uj b - 2u q 



A 2 A 2 



+ ui b - (2uj q - rj) ' 



(6) 
(7) 



which is very close to the exact value found by direct di- 
agonalization of the Hamiltonian (Fig. 2), and the effect 
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FIG. 2: (Color online) The frequency fizz for a truncated 
memory-qubit-bus system [the idling error is (flzzt) 2 ] for two 
values of the coupling: g m /2-K = gb/2-K = 25 MHz (blue lines) 
and g m /2n = gb/2n = 50 MHz (red lines). The solid lines 
show the results of exact diagonalization of the RWA Hamil- 
tonain (3) with and without g&. The effect of gs. is not visible 
(smaller than the line thickness). The blue and red dashed 
lines show the analytical result (6). 



of gd is of a higher order and therefore very small (see 
Appendix and Fig. 2). Notice that fizz <x V because in 
a linear system fizz = 0, and nonlinearity comes from 
the qubit. Equation (6) shows that an optimal choice of 
the qubit "parked" frequency is uj q = (u> m +u> b )/2, mid- 
way between the memory and the bus frequencies; then 
the idling error in this order goes to zero (this happens 
because the contribution of |020) in 1 101) becomes zero 
- see Appendix). Notice that in the RezQu architecture 
the frequencies of the memory resonators are assumed 
to be relatively close to each other (forming a "memory 
band" of frequencies). Then the optimal "parked" fre- 
quencies of the qubits are also close to each other. This 
is not a problem when all qubits are in state |0); how- 
ever, when a qubit is excited this may lead to a signif- 
icant resonant coupling with another qubit via the bus. 
To avoid this "spectral crowding" effect, it is useful to re- 
serve two additional frequencies, situated sufficiently far 
from the "parked" frequencies, at which a pair of qubits 
may undergo local rotations (simultaneous rotations of 
two qubits are often useful before and after two-qubit 
gates) . 

The idling error (2) scales quadratically with time. 
This is because we use a definition for which not the 
error itself but its square root corresponds to a metric, 
and therefore for a composition of quantum gates in the 
worst-case scenario we should sum square roots of the 
errors [37]. For the same reason, the worst-case idling 
error scales quadratically, Err cx A 2 , with the number 
N of qubits in a RezQu device. In principle, an average 
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idling error may scale linearly with N and time (for that 
we would need to define the memory "clock" frequency 
using an average occupation of the bus); however, here 
we use only the worst-case analysis. 

It is convenient to replace the time-dependence in the 
idling error estimate (NCtzz t) 2 by the dependence on 
the number of operations AT op in an algorithm. (The 
corresponding quadratic dependence on A^p can also be 
interpreted as the worst-case-scenario error for a compo- 
sition of quantum operations.) Assuming that each op- 
eration crudely takes time i op ~ g^ 1 + g^ 1 (this estimate 
comes from MOVE operations discussed later and also 
from controlled- Z gate) and neglecting the second factor 
in Eq. (6) (i.e. assuming non-optimal "parked" qubit fre- 
quencies) , we obtain the following estimate for the worst- 
case idling error: 



we obtain an error estimate Err ~ (N g^g 2 / A 4 ) 2 per op- 
eration. Since the worst-case scaling with the number of 
operations N op is always oc N% , we obtain the worst-case 
estimate 



Err : 



9tn9 b T N N U max( gm ,g b ) 
vam(g m ,g b ) ' 



A 4 A 4 



(8) 



where A m and Af> are typical dctunings at idling. Using 
for an estimate g m /27r = gb/^n = 25 MHz, i]/2ir = 200 
MHz, and A to /2tt = A fc /27r = 500 MHz, we obtain Err ~ 
lO^A^A^. 

To demonstrate the advantage of the RezQu architec- 
ture, Eq. (8) may be compared with the corresponding 
result for the conventional bus-based architecture (with- 
out additional memories). Then the idling error is due 
to ZZ-interaction between the qubit and the bus: the 
frequency of an idling qubit is affected by the bus occu- 
pation due to logic operations between other qubits. In 
this case ^zz.conv = -2g 2 ri/[A b (A b - 77)], that gives 



Err r . 



gWN 2 N 2 /A 



(9) 



Assuming a typical ratio g/ A < 0.1 between the coupling 
and detuning, we have a reduction in the idling error in 
the RezQu architecture by at least 10 even before con- 
sidering that Q,zz can be zeroed in this order. Using Eq. 
(9) we see that a conventional architecture allows only a 
very modest number of qubits and operations before the 
idling error becomes significant. In principle, the prob- 
lem can be solved by a constantly running dynamical 
decoupling (which would be quite nontrivial in a multi- 
qubit device). The RezQu idea eliminates the need for 
such dynamical decoupling. 

All our estimates so far were for the idling error due 
to the memory-bus interaction. Now let us discuss er- 
rors due to the four-step memory-memory interaction in 
the RezQu architecture. The ATAT-interaction between 
the memory and another (fcth) memory can be calcu- 
lated [39] as il X x = 2g m g m . k g b g b ,k/[A m A m . k (uj 
where additional subscript k indicates parameters for the 
fcth section of the device. The W-interaction does not 
produce a phase error accumulating in idling, but leads to 
the error Err ~ (Qxx / &m) 2 every time the information 
is retrieved from memory, where 5 m is the typical spacing 
between memory frequencies. Assuming similar sections 
of the RezQu device with ~ A r „ = A and S m ~ A/AT, 



Err- I^N^g^/A*. 



(10) 



This error is smaller than the idling error (8) if g m j, < V- 
For g m = gb = A/20 we find a very small error Err ~ 
10~ 10 N 2 N%, which means that there is essentially no 
spectral crowding problem for memories. Notice that for 
a conventional bus-based architecture the error estimate 
(10) is replaced by Err conv ~ N 2 N 2 p g 4 /A 4 and presents 
a difficult scaling problem due to the spectral crowding. 

Besides the XX-interaction between the two memo- 
ries, there is also the ZZ-interaction. Using the same 
approximate derivation as in Appendix [see Eq. (A7)] 
and assuming ~ A TO = A, we find an estimate 
Ozz « -r,g 4 m g 4 /A\ Then using Err ~ (n zz t op ) 2 N 2 p N 4 
(the scaling N 4 is because each pair brings a contribu- 
tion), with t op ~ g^ 1 + g^ , we obtain the worst-case 
estimate 



Err : 



jL9W n4n q P max(g m ,. 9fc ) 
A 16 mm(g m ,g b )' 



(11) 



This error is smaller than the memory-bus idling error (8) 
if N < A 4 /g^ n g%, which is always the case in practice. 



IV. MOVE OPERATION 

Any logic gate in the RezQu architecture requires mov- 
ing quantum information from one system element (mem- 
ory, qubit, bus) to another. Therefore the MOVE opera- 
tion is the most frequent one. It is important to mention 
that the one-way MOVE operation [41] is easier to design 
than the SWAP (or iSWAP) operation because we are 
not interested in the fidelity of the reverse transfer and 
can also assume zero occupation of the neighboring ele- 
ment. For example, for a perfect qubit— ^memory MOVE 
(zMOVE) operation in the truncated mqb system we 
search for a unitary, which transforms |010) — > — ij 100) 
(notice that in RWA |000) -> |000) always), but we are 
not interested in what happens to the initial states |001) 
and 1 100). Moreover, we can allow for an arbitrary phase, 
1 010) — > — ze~ tv |100), because this phase can be compen- 
sated either by shifting the operation start time within 
one period of the initial memory-qubit detuning [40] or 
by "qubit frequency excursion" with proper integral [7]. 
Therefore, we need to satisfy only two (complex) equa- 
tions to design the unitary [/move for this MOVE, 

(010|[/move|0T0) = 0, (00T|C/move|0T0) - 0. (12) 

Experimentally the qubit^memory MOVE is done 
[7, 15, 16, 40] by tuning the qubit in resonance (with some 
overshoot) with the memory resonator approximately for 
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a duration n/2g m . Equation (12) means that any rea- 
sonable shape of w q (t) tune/detune pulse with four ad- 
justable parameters can be used for a perfect MOVE op- 
eration in the truncated mqb system. Actually, as will 
be discussed later, the use of only two adjustable param- 
eters is sufficient to obtain an exponentially small error 
in the quasi-adiabatic regime. Such two-parameter con- 
struction is most convenient for practical purposes, but 
formally it is imperfect (non-zero error). So we will first 
discuss the perfect (zero error) four-parameter construc- 
tion. 

We have designed the qubit— » memory MOVE pulses 
Lo q (t) for the truncated mqb device both analytically (in 
first order) and numerically. The initial and final fre- 
quencies of the qubit are allowed to be different. In the 
analytical design we do calculations in the bare basis, 
|^(t)) = a(t) 1 100) +/3(t)|010) +7(t)|001), but define the 
co-moving frame as 



a e 



ft = e * IS fj = 7 e iUbt . (13) 



In this representation the only interesting initial state 
1 010) of the qubit— > memory MOVE is (in first order) 

5(0) = -WA m (0), m = 1, 7(0) = VA 6 (0), (14) 



where A r 



-ui q and Aft 



-u>b are the detunings. 



The desired (target) final state at time if is 



i.e. 



a(tf) 
0(*f) 



-Itp 



-lip 



, 7(if)=0, 



-i/ tf A m (t')dt' 



'| 100), 



(15) 



Notice that even though the phase tp is arbitrary, the 
relative phase between S(if) and /3(if) is fixed by the 
absence of the relative phase between a (if) and ft(tf). We 
see that the MOVE operation should eliminate the initial 
"tail" 7(0) on the bus (this needs two real parameters in 
the pulse design) and transfer most of the excitation to 
the memory with correct magnitude and relative phase 
of /3(if) (two more real parameters). 

Similarly to the experimental pulse design [7, 15, 16, 
40], we assume that the shape of uJ q [t) pulse consists of a 
front ramp, rear ramp, and a flat part in between them 
(Fig. 3 illustrates a piecewise-linear construction of the 
pulse) . As will be shown below, using two parameters for 
the front ramp shape we can ensure elimination of the 
"tail" 7(if); we can choose a rather arbitrary rear ramp, 
and using two parameters for the flat part (its frequency 
overshoot and duration) we can provide proper /3(if). 

Let us start with the "tail" 7. As follows from the 
Schrodinger equation with the Hamiltonian (3), 

7(i f ) = 7(0) - igt C dt~p{t)e- lS ° Ab(t ' )dt ' ■ (16) 



Let us denote the end of the front ramp by t\ and the 
start of the rear ramp by ti (see Fig. 4). For < t < t\ 
in Eq. (16) we can replace ft(t) with /3(0) = 1 because 
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FIG. 3: (Color online) Illustration of a piecewise-linear 
tune/detune pulse shape (qubit frequency as a function of 
time) for the MOVE operation qubit— ^memory in a three- 
component mqb system. The front ramp consists of two 
straight segments. The solid blue line shows the result of 
a four-parameter numerical optimization, in which the slope 
of the first straight segment, the qubit frequency at the end 
of the first straight segment, the duration — ii of the flat 
part and its overshoot D have been optimized. This gives 
Err = up to machine accuracy. The red dashed line shows 
analytical design based on Eqs. (17) and (22); in this case, 
Err = 5 x 10~ 4 . System parameters: tj m /27r = 7 GHz, 
uj b /2n = 6 GHz, u) q (Q)/2-K = 6.7 GHz, w g f»/27r = 6.5 GHz, 
g m /2ir — gt/2n — 25 MHz, the slopes of the second front seg- 
ment and of the final ramp have been fixed at 500 MHz/ns. 



the qubit occupation cannot change much during a short 
ramp. For t\ < t < i2 we can use integration by parts 
using Afc(t) fa u> m —u}f,i with ft(t) changing approximately 
from 1 to 0. Finally, there is a negligible (second-order 
in gb) contribution to the integral for t% < t < if because 
ft(t) is already small (first order in <?b). Thus for the 
desired pulse shape, in first order we obtain 



= 



7(*f) 
9b 



1 



A b (0) Vo ' 



-*Adt - 



(17) 



t" 



J% = I A b {t)dt. 



(18) 



As we see, required elimination of the "tail" 7 on the bus 
gives two equations (real and imaginary parts) for the 
front ramp shape. This can be done by using practically 
any shape with two adjustable parameters. Notice that 
condition (17) essentially means that in order to have 
correct (zero) "tail" on the bus at final time if, this tail 
at time if should be the same as the tail for the co-moving 
eigenstate of the qubit-bus system. 

Now let us design the flat part of the pulse, which 
should give us the proper ratio (3(tf)/a(t t[ ) from Eq. (15). 
After designing the front ramp we know a and ft at the 
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start of the flat part t\. in first order, P(t\ 

>S(*l) " A m (0) 



1 and 



e lB °dt, < 



t" 



A m (t)dt. 



Similarly, for an arbitrarily chosen rear ramp shape we 
know desired a and (3 at the end of the flat part t 2 : in 
first order, a(t 2 ) = 5(tf) = —ie~ lv and 



gfa) = 

a(t 2 ) A ro (t f ) 



e °° +ig m e 



e lB 'dt. 



(20) 



During the flat part of the pulse we can use the two- 
level approximation with coupling g m , and essentially 
connect the two points on the Bloch sphere correspond- 
ing to Eqs. (19) and (20) by a "Rabi" pulse. These points 
are close to the North and South poles, so the pulse is 
close to the ideal 7r-pulse; we assume a small constant 
overshoot A m = —D with \D/g m \ -C 1 (Fig. 4), and 
duration £ 2 — ^i = ^l^n ~ T with |r| <C "k/wr, where 
lur = \J 4g^ + D 2 . Then using the leading-order rela- 
tion for an almost perfect 7r-pulsc, 



p(t 2 ) a(h) D a(t 2 ) 

ig m T, ^— = -le 



rD/4g„ 



&{t 2 ) (3(h) 2g m ■ ~" ■ p(ti) 
we obtain the needed pulse parameters D and r as 



D 1 e~ lB o 

1- ir = 1 

2gl + A m (0) + A m (i f ) 

ti 



(21) 



+ i I e iB °dt + ie- lB o / e iB 'dt, (22) 



and also find the resulting phase ip = itD/Ag m . 

We have checked numerically the analytical pulse de- 
sign given by Eqs. (17) and (22). For example, for a 
piecewise-linear pulse whose front ramp consists of two 
straight segments (Fig. 4), the error 



Err=l-|(100|[/ M ovE|010)| ; 



(23) 



for the analytically designed pulses is found to be below 
10~ 3 for typical parameters with <ji m /27r = gb/^TT = 25 
MHz. As expected, the numerical four-parameter op- 
timization of such pulse shape gives zero error (up to 
machine accuracy) , and the shape of this perfect pulse is 
close to the analytically-designed shape (see Fig. 4). 

Experimental pulses for the MOVE operation [7, 15, 
16] are produced by a Gaussian filter and therefore have 
the error-function-shape ramps. We can use the same 
design idea for such pulses: shaping the front ramp us- 
ing two parameters (see Fig. 4) takes care of the "tail" 
7 (if) on the bus, for the rear ramp we use any conve- 
nient shape, and for the middle part we vary the over- 
shoot frequency and duration to ensure proper popula- 
tion transfer between the qubit and the memory (for such 
pulses it is natural to define the duration to be between 
the inflection points of the error- function shapes). We 
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FIG. 4: (Color online) Implementation of the qubit— ^memory 
MOVE operation in a three-component mqb system using a 
pulse with error-function-shaped ramps (sum of two time- 
shifted error functions for the front ramp). Four parameters 
of the pulse shape (see upper panel) are optimized: the time 
shift and amplitude ratio for the front-ramp error functions, 
the duration of the middle part of the pulse, and the overshoot 
magnitude. The error functions are produced by integrating 
Gaussians with standard deviation a = 1 ns; the beginning 
and end of the pulse are at 3u from the nearest error-function 
centers (shown by vertical lines in the upper panel). The 
middle and lower panels show time-dependence of the level 
populations in the bare-state basis and comoving eigenbasis. 
The MOVE error (23) is zero up to machine accuracy. 



have checked numerically that these four parameters are 
sufficient to achieve zero error (perfect transfer fidelity 
F = 1 — Err) in the truncated three-element system. 

As a further simplification of the MOVE pulse design, 
let us optimize only two middle-part parameters (over- 
shoot and duration) and do not optimize the front ramp 
shape. In this case we cannot ensure the proper "tail" 
7(tf); however, it is small by itself, and therefore the er- 
ror is not large. Moreover, for sufficiently slow pulses the 
"tail" is almost correct automatically because of the adi- 
abatic theorem. It is important to notice that the bus 
is well-detuned, \Ab/gb\ ^> 1, and then the adiabaticity 
condition is \dAj,/dt\ <C |A|/<7(,|, which is well-satisfied 
even by rather fast pulses. To estimate the correspond- 
ing error, we consider the two-level system bus-memory 
during the front ramp and write the differential equation 
for the variable y = 7//? — g^/A^,, which describes devia- 
tion from the co-moving eigenstate. Assuming \ j//3\ <C 1 
and |w g 7//3 2 | <C 1, we obtain y = iA^y + gbA b ~ 2 uj q . The 
"tail" error at the end of the front ramp is |y(ii) 2 | (no- 
tice that \(3(ti ) I w 1), and it does not change significantly 
during the rest of the pulse. If this is the major contri- 
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bution to the MOVE error, then 



Err = 



9b 



di q (t) exp(- 



-iAl) 



dt 



(24) 



where Aq is defined in Eq. (18). Numerical optimization 
of only the middle part of the pulse (overshoot and dura- 
tion) confirms that Eq. (24) is a good approximation for 
the MOVE error in this case. Notice that for an error- 
function ramp obtained by integrating a Gaussian with 
the standard deviation (time-width) a, the error (24) de- 
creases exponentially with a (we assume sufficiently long 
ramp time t\) and is typically quite small. For exam- 
ple, for Afc/27r changing from 0.5 GHz to 1 GHz and 
gb/2n = 50 MHz, the error is below 10~ 4 for a > 0.5 ns 
(for a > 0.35 ns if g h /2TT = 25 MHz). 

So far we have only considered the MOVE 
qubit— >-memory. The MOVE in the opposite direc- 
tion memory^qubit can be designed by using the time- 
reversed pulse shape. Perfect MOVE still requires opti- 
mizing four parameters (overshoot and duration of the 
middle part and the two parameters for the rear ramp), 
while using only the two parameters for the middle part 
is sufficient for a high-fidelity MOVE. 

In designing the MOVEs between the qubit and the 
memory we assumed no quantum information on the bus. 
The presence of an excitation on the bus makes the pre- 
viously designed perfect MOVE imperfect. We checked 
numerically that the corresponding error for typical pa- 
rameters is about 10 -4 , i.e. quite small. Moreover, a 
typical RezQu algorithm never needs a MOVE between a 
qubit and memory with occupied bus, so the unaccounted 
MOVE errors due to truncation are even much smaller. 
The analysis in this paper assumes the RWA Hamiltonian 
(3), which neglects terms cx a^jcr* and oc a^er^, which 
change the number of excitations. We have checked nu- 
merically that addition of these terms into the Hamil- 
tonian leads to negligibly small changes of the system 
dynamics during the MOVE operations. 

Designing MOVEs between the qubit and the bus is 
similar to designing MOVEs between the qubit and mem- 
ory, if we consider the truncated three-element system. 
However, in reality the situation is more complicated be- 
cause the bus is coupled with other qubits. Our four- 
parameter argument in this case does not work, and de- 
signing a perfect single-excitation MOVE would require 
2iV+2 parameters (for a truncated system with N qubits, 
one memory, and the bus), which is impractical. How- 
ever, the occupation of additional qubits is essentially the 
effect of the "tails" (if the discussed below problem of 
level crossing is avoided). Therefore, the desired "tails" 
can be obtained automatically by using sufficiently adi- 
abatic ramps in the same way as discussed above (for 
a MOVE qubit^bus the front ramp will be important 
for the "tails" from both sides, i.e. on the memory and 
other qubits). In analyzing the dynamics of the "tails" at 
other qubits, it is useful to think it terms of cigenstates 
of a truncated system, which includes the bus and other 
qubits (while excluding the qubit involved in the MOVE) . 



Then the "tail" error is the occupation of the eigenstates, 
mainly localized on other qubits. Since the frequencies 
of the bus and other qubits do not change with time, for 
the error calculation it is still possible to use Eq. (24) , in 
which Ab is replaced with tu q —uj q ,k (for the "tail" on fcth 
qubit), gb is replaced with gbg b ,k/(^ q ,k — and Aq is 
replaced with L [ui q (t') — uj q .k] dt' , where the subscript k 
labels additional qubit. This gives the estimate 



Err 



9bgb,kd) q {t) 
^lS)^b,k 



-ifo A,,*(t')dt' ^ 



(25) 



of the error due to the "tail" on fcth qubit for the 
qubit— >bus MOVE, in which integration is within the 
front ramp, A bjfc = u q ^ k -w b , and A qtk (t) = ui q {t) - w g;fe . 
The formula for the bus^qubit MOVE is similar, but 
the integration should be within the rear ramp. The er- 
ror (25) should be summed over additional N — 1 qubits 
(index fc) and therefore can be significantly larger than 
in our calculations for a truncated mqb system; however, 
the error increase is partially compensated by smaller ef- 
fective coupling g b gb,k/ ^-b.k- Crudely, we expect errors 
below 10 -4 for smooth ramps of few-nanosecond dura- 
tion and N < 10 2 . We emphasize that this simple solu- 
tion of the "tail" problem is possible only when we use 
eigenstates to represent the logical states. 

Another problem which we did not encounter in the 
analysis of the truncated system is the level crossing 
with other (empty) qubits during the MOVE operation. 
A simple estimate of the corresponding error is the fol- 
lowing. Effective resonant coupling between the moving 
qubit and another (fcth qubit) via the bus is gb9b,kl 'Ab, 
where gb and gb.k are two qubit-bus couplings and the 
detuning A b is the same for both qubits at the moment 
of level crossing. Then using the Landau-Zener formula 
we can estimate the error (population of the other qubit 
after crossing) as Err ~ 2Kg\g\ fe /(A 2 |w g |), where \io q \ 
is the rate of the qubit frequency change at the cross- 
ing. Our numerical calculations show that this estimate 
works well, though up to a factor of about 2 [when cur- 
vature of oj q {t) at the point of crossing is significant]. 
Using this estimate for g b /2n — g b ,k/^ — 25 MHz, 
A b /27r = 500 MHz, and uj q /2n = 500 MHz/ns, we ob- 
tain a quite significant error of about 10 -4 . A possible 
way to compensate this error is by using interference of 
the Landau-Zener transitions [42] on the qubit return 
transition. Another solution of the problem is to park 
empty qubits outside the frequency range between the 
bus and memories (above 7 GHz in our example). This 
would make impossible to cancel the idling error of Eq. 
(6) by using the "midway parking", but the idling er- 
ror is still small even without this cancellation [see the 
estimate below Eq. (8)]. Besides the qubit-qubit level 
crossings, there are also level crossings between a mov- 
ing qubit and other memories. This is a higher-order 
(weaker) process because of three steps between the qubit 
and memory. The effective coupling with fcth memory is 
then 3636,fcg m ,fc/(A b A mifc ), and the level crossing error 
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estimate is Err ~ 2irg 2 gl k g 2 nk / (A 2 A 2 mk \uj q \). 

In this paper we do not analyze two-qubit gates. Our 
preliminary numerical simulation of the controlled- Z gate 
has shown possibility of a high-fidelity gate design (with 
the error of about 10~ 3 , mainly due to level crossing). 
However, we have not studied this gate in detail. A de- 
tailed analysis of two-qubit gates in the RezQu architec- 
ture will be presented elsewhere [36] . 



V. TUNNELING MEASUREMENT 

Finally, let us discuss whether or not using the eigen- 
states as the logical states presents a problem for mea- 
surement. Naively, one may think about a projective 
measurement of an individual qubit; in this case the logic 
state "1" would be erroneously measured as "0" with 
probability of about (g m /A m ) 2 + (gb/Ab) 2 ~ 10~ 2 be- 
cause the eigenstate spreads to the neighboring memory 
and bus. This would be a very significant error, and 
the bare-state representation of logical states would be 
advantageous. However, this is not actually the case be- 
cause any realistic measurement is not instantaneous (not 
projective). In fact, if a measurement takes longer than 



(e.g., [43]) 



then the eigenbasis is better than the bare basis. 
As a particular example let us analyze tunneling mea- 
surement of a phase qubit [14] (we expect a similar re- 
sult for the qubit measurement in the circuit-quantum- 
electrodynamics (cQED) setup [26, 31]). The bare states 
|0) and |1) of a phase qubit correspond to the two lowest 
energy states in a quantum well, and the measurement 
is performed by lowering the barrier separating the well 
from essentially a continuum of states [14]. Then the 
state |1) tunnels into the continuum with a significant 
rate T, while the tunneling rate for the state |0) is negli- 
gible. The event of tunneling is registered by a detector 
"click" (the detector is a SQUID, which senses the change 
of magnetic flux produced by the tunneling) . In the ideal 
case after waiting for a time t » the measurement 
error is negligibly small (in real experiments the ratio of 
the two tunneling rates is only ~ 10 2 , which produces a 
few-per-cent error; however, we neglect this error because 
in principle it can be decreased by transferring the state 
|1) population to a higher level before the tunneling, and 
also because here we are focusing on the effect of "tails" 
in the neighboring elements). 

In presence of the memory and resonator coupled to 
the qubit, the logic state "0" still cannot be misidenti- 
fied, because the tunneling is impossible without an exci- 
tation. However, the logic state "1" can be misidentified 
as "0" , when sometimes the expected tunneling does not 
happen (because part of the excitation is located in the 
memory and resonator). Let us find probability of this 
error. For simplicity we consider a two-component model 
in which a phase qubit is coupled to its memory resonator 
only, and restrict the state space to the single-excitation 
subspace of this mq system. Then the tunneling pro- 
cess can be described by the non-Hermitian Hamiltonian 



H = 



iT/2 



(26) 



and the error in measuring the logic state "1" (identifying 
it as "0" after measurement for time t) is its survival 
probability 



Ew=\m)\m)\ 2 , 



(27) 



where the initial state is normalized, | (-0(0) 1-0(0)) | 2 = 1- 
Our goal is to compare this error for the cases when the 
initial state is the bare state 10(0)) = 1 01) or the eigen- 
state |0(O)) = 1 01) (in this mq notation the qubit state 
is shown at the second place, and 1 01) is the eigenstate 
before the measurement, i.e. when T = 0). 

The solution |0(i)) of the time-dependent Schrodingcr 
equation, i(d/dt)\tp(t)) = H\ij)(t)), is given by the lin- 
ear combination, \^{t)) = C m \W) e - lE ^ + C q \01) e"^*. 
Here the eigenstate notation with the tilde sign reminds 
of a non-zero T, the constants C m ^ q depend on the ini- 
tial conditions, and E„^ q = Re(E m ^ q ) — iY mA j2 are the 
complex eigenenergies, which include the corresponding 
decay rates r m and T q of the eigenstates located mainly 
on the memory and the qubit. Diagonalizing the Hamil- 
tonian (26) and assuming weak coupling, g m -C A m , 
g m < T, we find 



A 2 „ + (r/2) 2 



«r, 



r — r 



r. (28) 



For measurement during a sufficiently long time i > T, 
only the |10)-term in \ip(t)) survives, and correspond- 
ingly the error (27) is Err = |C„ 
(Toj0(O)>. Thus we obtain 



|2 e r m t w h ere n 



Err P 



Err 



bare 



4Ai 



Err 



bare 



(29) 



for the measurement errors starting either with the eigen- 
state or with the bare state. Even though both errors 
decrease with the measurement time t as e~ Tmt , the rate 
r r „ is small [see Eq. (28)], so for a realistically long mea- 
surement we can use e~ Tmt w 1. 

Equation (29) shows that from the measurement point 
of view it is advantageous to use the eigenstates to rep- 
resent the logical states rather than the bare states if 
r < 2A m . For a typical value A m /27r = 0.5 GHz this 
requires > 0.16 ns, which is always the case. 

Figure 4 shows the state dynamics during the tunnel- 
ing measurement in the bare basis \ip{t)) = a(t) 1 10) + 
j3(t) 1 01) , starting either with the eigenstate 1 01) or with 
the bare state 1 01) . The oscillations correspond to the 
beating frequency A m /27r = 0.5 GHz. We see that sim- 
ilarly to the above-analyzed dynamics in the eigenba- 
sis, |/3(i)| 2 becomes exponentially small after t 3> 
while |a(i)| 2 essentially saturates (decaying with a much 
smaller rate T m ). Then for the assumed tunneling 
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FIG. 5: (Color online) Time dependence of the squared am- 
plitudes |a| 2 and \f3\' 2 of the mq state \ip(t)) = a(t)\10) + 
fi(t) 1 01} , decaying in the process of tunneling measurement. 
Blue curves correspond to the system initially prepared in the 
bare state IV'(O)} = |01), while for red curves the initial state 
is the eigenstate | i/j (0) } = 1 01) . For t 2> T _1 the measure- 
ment error (27) is mainly the residual occupation |o;(t) | 2 of 
the memory resonator (solid curves). For the depicted system 
parameters |a e igen(£)| 2 /|«bare(£)| 2 ~ 0.025 at t > 10 ns, i.e. 
the error of the eigenstate measurement is 40 times smaller 
than the error of the bare state measurement. 



the RezQu devices seems to be needed. 

We have shown that instead of using bare states it 
is much better to use eigenstates to represent logical 
states. In this case there is essentially no dynamics in 
idling (except for the phase errors), which greatly sim- 
plifies a modular construction of a quantum algorithm. 
The logical encoding by eigenstates is also advantageous 
for the single-qubit state generation and measurement. 
We have presented a simple design for the MOVE opera- 
tion, which is the most frequent operation in the RezQu 
architecture. We have shown that a four-parameter op- 
timization is sufficient for designing a perfect MOVE in 
a truncated three- component system. Moreover, opti- 
mization of only two experimentally-obvious parameters 
is sufficient for high-fidelity MOVEs (with errors less than 
10~ 4 ). While we have not analyzed two-qubit gates, we 
expect that their design with similar high fidelity is also 
possible. Overall, we believe that the RezQu architecture 
offers a very significant advantage compared to the previ- 
ously proposed architectures for superconducting qubits, 
and we believe that this is the practical way to progress 
towards a medium-scale quantum computing device. 

This work was supported by IARPA under ARO grant 
W911NF-10-1-0334. The authors thank Michael Geller, 
Farid Khalili, Matteo Mariantoni, Leonid Pryadko, and 
Frank Wilhelm for useful discussions. 



Appendix A: DERIVATION OF Vlzz 



1 the ratio of errors Err G i gon /Errb a 

>are 



rate T = 1 ns~ 

IcfceigenWlVlckbaret^)! 2 (with the subscript denoting the 
initial state) saturates at approximately the value 0.025 
given by Eq. (29) . 

We emphasize that even though we have only consid- 
ered the tunneling measurement, the result (29) for the 
measurement error is expected to remain crudely valid 
for most of realistic (i.e. "weak") measurements with a 
time scale In particular, for the cQED setup we 

expect that the role of T is played (up to a factor) by the 
ensemble dephasing rate due to measurement. 



In this Appendix we derive Eq. (6) for fizz = £101 + 
£ ooo~ eioo — eooi in the truncated mqb system. We assume 
that the couplings g m and gb are of the same order, g m ~ 
gb ~ g, and do calculations in fourth order in g. 

The RWA Hamiltonian (3) leads to the forma- 
tion of three subspaces, which do not interact with 
each other: the ground state |000) = 1 000) (with 
zero energy, eooo = 0), the single-excitation subspace 
{ 1 100) , |010), |001)}, and the two-excitation subspace 
{|101), jllO), joil), 1 200), |020), |002)}. It is rather easy 
to find eigenenergies in the single-excitation subspace; 
neglecting the direct coupling gd, in fourth order in g we 
obtain 



VI. CONCLUSION 

In summary, we have discussed the main ideas of the 
RezQu architecture and analyzed several error mecha- 
nisms, excluding analysis of two-qubit gates. The main 
advantage of the RezQu architecture is the strong (> 10 4 
times) reduction of the idling error compared to the con- 
ventional bus-based architecture, and also an effective so- 
lution to most of the problems related to spectral crowd- 
ing. In the absence of decoherence this makes possible a 
simple scaling of a RezQu device to ~ 30 qubits without 
the need for dynamical decoupling. For further scaling 
the next architectural level of communication between 



e ioo 



eooi — 



Srr. 



9m 

A 2 



2 2 
9m9b 



9 m 9b 



Ag(w m -u>b) 



(Al) 
(A2) 



To find eioi, we write the eigenstate 1 101) as a super- 
position of all elements of the two-excitation subspace, 

|T01) = [|101> +a U o|H0) +ooii|011> +a 200 |200) 

+a 020 |020) + a o 2 |002)] /Norm (A3) 

with unimportant normalization. Then the Schrodinger 
equation iJ|101) = eioi 1 101) (again neglecting gj) gives 
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six equations: 

(u m + L0 b ) + a lw g b + a 01 ig m = eioi, 

9b + aiio(w m + w 9 ) + (o 2 oo + a 020 )g m V2 = on eioi, 

9m + aou(u q +uj b ) + (a 020 + a 002 )g b V2 = a iieioi, 

m ^200^101; 

anog m V2 + a Q ng b V2 + a Q2Q (2uj q - rf) = a 2oeioi, 
aoiiffbV2 + a 002 2uj b = a Q02 e wl . (A4) 

From the first three of them we obtain 

9 m + ( a 020 + a 002 )g m g b y/2 



£101 ~ 



+ 



1 + Wfe ' 

9 b + (^200 + Q02o)gmgbV / 2 
£101 — w m — w„ 



(A5) 



which gives eioi in fourth order in g if we use second- 
order eioi in the denominators (which is obtained from 
the same equation using zeroth-order eioi) and second- 
order amplitudes 0200, 0:020 > 0002 ■ These amplitudes can 
be found from the last three equations (A4) using the 
first-order values a 110 = -g b /A b , a 011 = g m /A m : 



g m gbV2 g m gb\/2 

O200 — "7 7 r, O002 — 



O020 — 



A b (u; m -u; b )' A m (w m -u> 6 ) 

-g m gbV2 uj m + u b - 2ui q 
A m A b uj m + u b - (2u q - 77) 



(A6) 



Finally, substituting Eq. (A6) into Eq. (A5), and using 
Eqs. (Al), (A2), we obtain Eq. (6) for fl Z z in fourth 
order. 

The above was the formal derivation of Eq. (6). Let 
us also obtain it approximately. Since in a linear system 
Ozz = (excitations do not interact with each other), 
a non-zero value can come only from the qubit nonlin- 
earity r\. Assuming small 77, we can use the first-order 



perturbation theory in 77 to find the energy shift of 1 101) 
due to contribution from 1 020) (occupation of the qubit 
second level), 



zz 



Se 10 i 



-V I "020 1 



(A7) 



To find «020 we start with the first-order (in g) eigenstate 
1 101) = |101) - (g b /A b ) |110) + (g m /A m ) |011) and then 
obtain the next-order estimate 0:020 = ^J2g m g b \ij ' A b — 
l/A m ]/(2(x> 9 — 77 — uj m — uj b ) [which coincides with Eq. 
(A6)]. If this estimate is substituted into Eq. (A7), then 
for Slzz we obtain Eq. (6) with the squared second frac- 
tion. However, if in the above formula for 0020 we neglect 
r] (as we should for first-order perturbation in 77), then 
we obtain Eq. (6) with the second fraction replaced by 
1. One may say that the average between these two re- 
sults for small 77 reproduces Eq. (6); however, it is more 
appropriate to say that this approximation (first order in 
77) can accurately reproduce only the first fraction in Eq. 
(6), while the second fraction is beyond the accuracy of 
the approximation. 

A slightly different derivation reproduces Eq. (6) ex- 
actly. Instead of using the first-order approximation in 
77 [Eq. (A7)], let us find the change of eioi due to its 
repulsion from e 2o- Since the effective interaction is 
<7off = V2g m g b [l/ A b — 1/A m ] (see the above estimate of 
0:020), the repulsion is 5e 10 i = -g^ s /(2uj q -r]-uj m -ui b ). 
The difference between this repulsion with and without 
nonlinearity 77 gives Q. Z z = -5eff[( 2w « -V-Um -^b) -1 - 
(2uj q — ui m — Wfc) -1 ], which reproduces Eq. (6). 

These simple derivations do not change if we take into 
account the direct interaction gd in the Hamiltonian (3), 



which is of the second order, g^ 



because of Eq. 



(5). Therefore, the fourth-order result (6) for VLzz should 
not change either. The rigorous fourth-order calculation 
shows that ei 00 + £001 increases by 2g d g m g b /A m A b due 
to gd, but eioi increases by the same amount, so that 
these contributions to Sl Z z cancel each other. 



[1] J. Clarke and F. K. Wilhelm, Nature 453, 1031 (2008). 
[2] T. Yamamoto, Yu. A. Pashkin, O. Astafiev, Y. Naka- 

mura, and J. S. Tsai, Nature 425, 941 (2003). 
[3] J. H. Plantenberg, P. C. de Groot, C. J. P. M. Harmans 

and J. E. Mooij, Nature 447, 836 (2007). 
[4] M. Steffen, M. Ansmann, R. C. Bialczak, N. Katz, E. 

Lucero, R. McDermott, M. Neeley, E. M. Weig, A. N. 

Cleland and J. M. Martinis, Science 313, 1423 (2006). 
[5] L. DiCarlo, J. M. Chow, J. M. Gambetta, L. S. Bishop, 

B. R. Johnson, D. I. Schuster, J. Majer, A. Blais, L. 
Frunzio, S. M. Girvin and R. J. Schoelkopf, Nature 460, 
240 (2009). 

[6] L. DiCarlo, M. D. Reed, L. Sun, B. R. Johnson, J. M. 

Chow, J. M. Gambetta, L. Frunzio, S. M. Girvin, M. H. 

Devoret and R. J. Schoelkopf, Nature 467, 574 (2010). 
[7] M. Mariantoni, H. Wang, T. Yamamoto, M. Neeley, R. 

C. Bialczak, Y. Chen, M. Lenander, Erik Lucero, A. D. 



O'Connell, D. Sank, M. Weides, J. Wenner, Y. Yin, J. 
Zhao, A. N. Korotkov, A. N. Cleland, J. M. Martinis, 
Science 334, 61 (2011). 
[8] T. Hime, P. A. Reichardt, B. L. T. Plourde, T. L. Robert- 
son, C. E. Wu, A. V. Ustinov, and J. Clarke, Science 314, 
1427 (2006). 

[9] F. Altomare, J. I Park, K. Cicak, M. A. Sillanpaa, M. S. 
Allman, D. Li, A. Sirois, J. A. Strong, J. D. Whittaker, 
and R. W. Simmonds, Nature Phys. 6, 777 (2010). 

[10] J. Bylander, S. Gustavsson, F. Yan, F. Yoshihara, K. 
Harrabi, G. Fitch, D. G. Cory, Y. Nakamura, J. S. Tsai, 
and W. D. Oliver, Nature Phys. 7, 565 (2011). 

[11] F. Mallet, F. R. Ong, A. Palacios-Laloy, F. Nguyen, P. 
Bertet, D. Vion, and D. Esteve, Nature Phys. 5, 791 
(2009). 

[12] R. Vijay, D. H. Slichter, and I. Siddiqi, Phys. Rev. Lett. 
106, 110502 (2011). 



11 



[13] For a high-level architecture solution, see D. P. DiVin- 

cenzo, Phys. Scr. T137, 014020 (2009). 
[14] J. M. Martinis, Quant. Inf. Processing 8, 81 (2009). 
[15] H. Wang, M. Mariantoni, R. C. Bialczak, M. Lenander, 

E. Lucero, M. Neeley, A. O'Connell, D. Sank, M. Wei- 

des, J. Wenner, T. Yamamoto, Y. Yin, J. Zhao, J. M. 

Martinis, A. N. Cleland, Phys. Rev. Lett. 106, 060401 

(2011). 

[16] M. Mariantoni, H. Wang, R. C. Bialczak, M. Lenander, 
E. Lucero, M. Neeley, A. D. O'Connell, D. Sank, M. Wei- 
des, J. Wenner, T. Yamamoto, Y. Yin, J. Zhao, J. M. 
Martinis, A. N. Cleland, Nature Phys. 7, 287 (2011). 

[17] S. T. Merkel and F. K. Wilhelm, New J. Phys. 12, 093036 
(2010). 

[18] F. W. Strauch, P. R. Johnson, A. J. Dragt, C. J. Lobb, 

J. R. Anderson, and F. C. Wellstood, Phys. Rev. Lett. 

91, 167005 (2003). 
[19] G. Haack, F. Helmer, M. Mariantoni, F. Marquardt, and 

E. Solano, Phys. Rev. B 82, 024514 (2010). 
[20] T. Yamamoto, M. Neeley, E. Lucero, R. C. Bialczak, J. 

Kelly, M. Lenander, M. Mariantoni, A. D. OConnell, D. 

Sank, H. Wang, M. Weides, J. Wenner, Y. Yin, A. N. 

Cleland, and J. M. Martinis, Phys. Rev. B 82, 184515 

(2010). 

[21] A. O. Niskanen, K. Harrabi, F. Yoshihara, Y. Nakamura, 

5. Lloyd, and J. S. Tsai, Science 316, 723 (2007). 

[22] R. C. Bialczak, M. Ansmann, M. Hofheinz, E. Lucero, M. 
Neeley, A. D. O'Connell, D. Sank, H. Wang, J. Wenner, 
M. Steffen, A. N. Cleland, J. M. Martinis, Nature Physics 

6, 409 (2010). 

[23] S.-L. Zhu, Z. D. Wang, and K. Yang, Phys. Rev. A 68, 
034303 (2003). 

[24] A. Blais, A. M. van den Brink, and A. M. Zagoskin, Phys. 

Rev. Lett. 90, 127901 (2003). 
[25] X. Zhou, M. Wulf, Z. Zhou, G. Guo, and M. J. Feldman, 

Phys. Rev. A 69, 030301 (2004). 
[26] A. Blais, R.-S. Huang, A. Wallraff, S. M. Girvin, and R. 

J. Schoelkopf, Phys. Rev. A 69, 062320 (2004). 
[27] A. N. Cleland and M. R. Geller, Phys. Rev. Lett. 93, 

070501 (2004). 

[28] A. Blais, J. Gambetta, A. Wallraff, D. I. Schuster, S. M. 
Girvin, M. H. Devoret, and R. J. Schoelkopf, Phys. Rev. 



A 75, 032329 (2007). 
[29] J. Koch, T. M. Yu, J. Gambetta, A. A. Houck, D. I. 

Schuster, J. Majer, A Blais, M. H. Devoret, S. M. Girvin, 

and R. J. Schoelkopf, Phys. Rev. A 76, 042319 (2007). 
[30] J. Majer, J. M. Chow, J. M. Gambetta, Jens Koch, B. 

R. Johnson, J. A. Schreier, L. Frunzio, D. I. Schuster, A. 

A. Houck, A. Wallraff, A. Blais, M. H. Devoret, S. M. 

Girvin and R. J. Schoelkopf, Nature 449, 443 (2007). 
[31] R. J. Schoelkopf, S. M. Girvin, Nature 451, 664 (2008). 
[32] J. Q. You, F. Nori, Nature, 474, 589 (2011). 
[33] E. J. Pritchett and M. R. Geller, Phys. Rev. A 72, 

010301(R) (2005). 
[34] M. A. Sillanpaa, J. I. Park, and R. W. Simmonds, Nature 

449, 438 (2007). 
[35] B. R. Johnson, M. D. Reed, A. A. Houck, D. I. Schuster, 

L. S. Bishop, E. Ginossar, J. M. Gambetta, L. DiCarlo, 

L. Frunzio, S. M. Girvin and R. J. Schoelkopf, Nature 

Physics 6, 663 (2010). 
[36] J. Ghosh et al., Controlled-Z logic gate in 

Resonator /zero-Qubit architecture, in preparation 

(2011). 

[37] M. A. Nielsen and I. L. Chuang, Quantum Computation 

and Quantum Information (Cambridge University Press, 

Cambridge, UK, 2000). 
[38] M. Boissonneault, J. M. Gambetta, and A. Blais, Phys. 

Rev. A 79, 013819 (2009). 
[39] R. A. Pinto, A. N. Korotkov, M. R. Geller, V. S. 

Shumeiko, and J. M. Martinis, Phys. Rev. B 82, 104522 

(2010). 

[40] M. Hofheinz, H. Wang, M. Ansmann, R. C. Bialczak, 
Erik Lucero, M. Neeley, A. D. O'Connell, D. Sank, J. 
Wenner, J. M. Martinis and A. N. Cleland, Nature 459, 
546 (2009). 

[41] L. P. Pryadko and A. N. Korotkov, unpublished. 

[42] W. D. Oliver, Yang Yu, J. C. Lec, K. K. Berggren, L. 

S. Levitov and T. P. Orlando, Science 310, 1653 (2005); 

M. Sillanpaa, T. Lehtinen, A. Paila, Yu. Makhlin, and P. 

Hakonen, Phys. Rev. Lett. 96, 187002 (2006). 
[43] K. Gottfried, V. F. Weisskopf, Concepts of Particle 

Physics, Vol. I (Oxford University Press, 1984), p. 151. 



